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ABSTRACT 

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux 
which governs their radiative and dynamical properties. We compute a suite of 3D circulation models 
with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with 
and without temperature inversions. We find that, for irradiation temperatures T- ln - ^, 2000 K, heat 
redistribution is very efficient, producing comparable day- and night-side fluxes. For Tj rr tzs 2200- 
2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our 
simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective 
to radiative timescales. Models with temperature inversions display a higher day-night contrast due 
to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary 
compared to the effects of irradiation. The hotspot offset from the substellar point is large when 
insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The 
atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the 
Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of 
a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic 
dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the 
penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases 
steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, 
favoring radius inflation in the most irradiated objects. A survey of the existing data, as well as the 
inferences made from them, reveals that our results are broadly consistent with the observational 
trends. 

Subject headings: planets and satellites: atmospheres 



1. INTRODUCTION 

Hot Jupiters are a class of exoplanets believed to have 
their rotational period synchronized with their orbital 
motion due to the strong tidal forces resulting from the 
proximity to their parent stars. Because of tidal locking, 
these exoplanets possess permanent day- and night-side 
hemispheres. These special conditions make hot Jupiters 
exciting astrophysical laboratories, in which the interplay 
of radiative forcing, atmospheric dynamics, composition, 
rotation speed and surface gravity contribute to their 
observational properties. 

A growing body of observations in recent years has 
allowed for hot Jovian atmospheres to be characterized 
via a combination of different techniques, which include 
the measurements of transmission spectra, thermal phase 
curves, and secondary eclipses. A complex picture has 
emerged with the objects characterized by a wide range 
of properties: some display temperature inversions in 
their atmospheres, others do not; heat redistribution is 
very efficient in some cases, but highly inefficient in oth- 
ers; a fraction of exoplanets display radii much too large 
to be compatible with standard evolutionary models. 

Theoretical explanations for this variety of phenomena 
have been abundant. The presence of temperature in- 
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versions in some exoplanets, but not in others, has been 
shown to be an effect of atmospheric opacity (Fortney 
et al. 2008). The variety of recirculation efficiencies, on 
the other hand, may be attributed to several mechanisms 
(e.g., Cowan et al. 2012), such as shorter radiative times 
in hotter exoplanets, stronger magnetic drag (which re- 
tards the wind speeds) at higher irradiation fluxes, and 
a weaker greenhouse with increasing exoplanetary tem- 
perature. 

For the problem of inflated radii, a variety of possible 
explanations has been put forward (see, e.g., Fortney & 
Nettelman 2010 for a recent review). Bodenheimer et al. 
(2001) pointed out that tidal dissipation, induced by a 
non-negligible orbital eccentricity, can pump energy into 
the interior and hence counteract contraction during evo- 
lution. This idea has been further explored and refined 
by Jackson et al. (2008), Miller et al. (2009) and Ibgui 
& Burrows (2009). Other ideas (specifically proposed 
for the case of HD 209458b) include the tidal damping 
of obliquity (Winn & Holman 2005) and extreme evapo- 
ration (Baraffe et al. 2004). Further, it was noted that 
enhanced opacities in the atmosphere may retard cool- 
ing and hence delay contract ion (Burrows et et al.. 2007a, 
but see also Guillot 2008). iGuillot fc Showman 1 ([20021 ) 
have suggested that a fraction of a percent or so of the 
absorbed stellar flux can be converted into kinetic energy 
in the deep interior by the breaking of atmospheric waves, 
but Burkert et al. (2005) were not able to confirm this 
dissipation mechanism via numerical simulations. More 
recently, the observation that the flow of ionized particles 
in the interior generates substantial atmospheric currents 
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has led to investigations which showed that Ohmic dissi- 
pation can inject enough energy in the deep atmospheric 
levels to affect the evolution of the radius (Batygin & 
Stevenson 2010; Perna et al 2010b). 

Ultimately, discrimination among various models will 
happen through a close comparison with the data, as 
more accumulate. From an observational standpoint 
alone, current data strongly hint at the fact that the 
strength of the irradiating flux plays a major role in de- 
termining the degree of heat redistribution (Harrington 
2011), as well as the amount of radius inflation (e.g., 
Fortney et al. 2010). This is unsurprising, since irra- 
diation is the main energy provider in these exoplanets. 
However, from a theoretical point of view, the direct link 
between strength of irradiation and those fundamental 
observational properties of the exoplanets has not been 
elucidated yet. This has motivated our current work. 

In this paper, by means of global, 3D atmospheric cir- 
culation models with dual-band radiative transfer ta i- 
lored to hot Jupiters (jHeng. Frierson fe Phillippj|2011| ). 
we explore, as a function of the irradiating flux from the 
parent star: 

a ) The efficiency of heat redistribution from day to night 
side, focusing on the infrared photosphere, and the rela- 
tive day-night brightness contrast; 

b) The offset of the hotspot (hottest infrared region of 
the exoplanet) with respect to the substellar point; 

c) The magnitude of hydrodynamic dissipation due to 
both the onset of the Kelvin-Helmholtz instability and 
shocks; 

d) The strength of Ohmic dissipation. 

We present a suite of 16 models: for eight values of 
the irradiation temperature (spanning the range T; rr ps 
770 — 3000 K), and each for cases with and without a 
temperature inversion. 

Our simulations allow us to elucidate the relative im- 
portance of irradiation and opacity (correlated with tem- 
perature inversions) in determining the observational 
properties of the exoplanet. They also allow us to assess 
what is the main effect in determining heat redistribu- 
tion and day-night contrast, and to predict theoretically, 
for the first time, the temperature or irradiating flux 
around which redistribution begins to break down. Fur- 
thermore, by computing the Richardson number and the 
Mach number at each point in the flow, we can quantita- 
tively estimate the depth down to which the flow can be 
subject to either shear instabilities or shocks, and hence 
whether hydrodynamic dissipation of this type can af- 
fect radius inflation. Similarly, our simulations allow us 
to follow the depth down to which Ohmic power can be 
dissipated, as a function of the irradiation strength, and 
for models with varying opacities. 

Our paper is organized as follows: in §2, we summa- 
rize the main features of our simulations, and provide 
the specific details for the models that we study. Our re- 
sults for heat redistribution and the hotspot location are 
presented in §3, while in §4 we study hydrodynamic and 
Ohmic dissipation. A qualitative comparison with obser- 
vations is made in §5. We finally summarize and conclude 
in §6. For completeness, we provide zonal-mean profiles 
of the zonal wind, temperature, potential temperature 
and streamfunction in Appendix [AJ 

2. METHODOLOGY 



Our 3D atmospheric circulation models 
are based on the computationa l setup de- 
scribed in iHeng. Menou fe P hillipps (|2011f) and 
iHeng. Frierson fc Phillippsl (|2011ft . In the follow- 
ing, we summarize its main features, as well as the set 
of parameters which we use for the current simulations. 

2.1. Stellar and orbital parameters 

A concise list of the parameters and their adopted val- 
ues are described in Table [1] Most of the parameters are 
chosen to have typical values, e.g., the surface gravity 
of the hot Jupiter is g p = 10 m s~ 2 . Since our primary 
interest is to explore the dependence of selected exoplan- 
etary properties on the strength of irradiation and the 
shortwave atmospheric opacity, we keep the rest of the 
parameter values to be the same for all of the models 
(Table P. 

Assuming a Bond al bedo of zero, the i rradiation^! tem- 
perature is defined as (|Heng et "ail 120 12ft 
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such that the irradiating flux is F — osbT£ t . (For 
parameter values appropriate to Earth, one recovers 
J~a ~ 1370 W m~ 2 , the solar constant.) The irradiat- 
ing flux J-"o, stellar radii i?*, effective stellar temperature 
T* and exoplanet-star spatial separation a are related by 
the expression, 



a = R±Tt 
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where <tsb is the Stefan-Boltzmann constant and we have 
assumed that the exoplanet orbits a Sun-like star. The 
use of the irradiation, rather than the equilibrium tem- 
perature, circumvents the issue of dealing with (confus- 
ing) factors of order unity related to assumptions about 

the efficiency of heat redistrib ution. 

The tidal locking time is (jGoldreich fc Sote"n 119661 : 
iBodenheimer et al. Il200lh 
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where uj p and Q p are the rotational and orbital frequen- 
cies of the exoplanet, respectively. For our models, we 
get t iock ~ 10 3 -10 8 yr (Q/lQ 6 ){w p /Q. p ). Thus, for a Gyr- 
old solar-like star (which is what we assume), our as- 
sumption of tidal locking is plausible and self-consistent 
if Wp/flp ~ 1. If uj p /fl p ~ 10, then this assumption be- 
comes suspect only for Model C. With the rotational and 
orbital frequencies being equal, we obtain 



(GM+ 
\ a 3 
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where G is the gravitational constant and M± is the stel- 
lar mass. Thus, J-o, Q p and a need to be varied self- 

5 In this study, we use the terms "irradiation" and "isolation" 
interchangeably. 



3 



TABLE 1 

Table of common parameters and their values 



Symbol 


Description 


Units 


Value(s) 




norrricilizBjtioii for optical depth of shortwave absorbers 




5 x 10 3 5 x 10 4 




normalization lor optical depth ot longwave absorbers 




-I r»4 

10 


e 


enhancement factor due to collision-induced absorption 




1 


cp 


specific heat capacity at constant pressure (of the atmosphere) 


J K- 1 kg" 1 


14550.4 


n 


specific gas constant (of the atmosphere) 


J K- 1 kg" 1 


4157.25 


K = 1Z/Cp 


adiabatic coefficient^ 


J K- 1 kg" 1 


2/7 


Po 


reference pressure at bottom of simulation domain 


bar 


1000 


9p 


acceleration due to gravity 


_2 

m s 


10 


Rp 


radius of hot Jupiter 


km 


7.1 x 10 4 


tu 


hyperviscous time 


dayt 


10" 6 


N v 


vertical resolution 




36 



t: expressed in terms of an exoplaneta ry day (i.e., rot ational period), 
we term n the "adiabatic coefficient", following Pierrchumbcrt (2010), and avoid calling it 
the "adiabatic index" in order to not confuse it with the ratio of specific heat capacities. 
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consistently. We note that setting the rotational and or- 
bital frequencies to be equal results in a more constrained 
problem for our simulations. While the existence of close- 
in, non-spin-synchronized gas giants remains possible, 
there is little empirical evidence currently available to 
support such a modelling effort. 

We present 16 models in this paper, divided into sub- 
groups which we term "cold" (C, CI and C2), "warm" 
(W, Wl and W2) and "hot" (H, HI). For Model W, 
the irradiating flux J-p is chosen to match the threshold 
value found by iDemorv fc Seagerl (|2011[ ). below which 
transiting Kepler giant exoplanet candidates appear to 
have non-inflated radii. Each of these models is stud- 
ied both with and without the presence of a temperature 
inversion. The parameters specific to each model are re- 
ported in Table H 

2.2. Thermodynamics 

The adiabatic coefficient k is related to the 
thermodynamical pr operties of the atmospheric gas 
(|Pierrehumbertll2010[ ). 

= ?*L - 2 

c P 2 + n dof 

By assuming the atmosphere to be dominated by molec- 
ular hydrogen (i.e., the number of degrees of freedom 
of the gas is n^of = 5), we get k — 2/7. Further- 
more, the specific gas constant is related to the mean 
molecular weight pL and the universal gas constant 1Z* by 
K = K*/n = 4157.25 J K" 1 kg- 1 . We use /i = 2 to 
emphasize the detachment from any specific case study, 
and remark that our results are rather insensitive to small 
changes in the mean molecular weight (e.g., /i = 2.35 as 
is typical for a solar mix). It follows that the specific heat 
at constant pressure is cp = 14550.4 J K _1 kg -1 . We 
recall that the adiabatic index is 7 = 1 + 2/ridof = 7/5 
for an atmosphere dominated by molecular hydrogen and 
that the sound speed is c s = \J r ykBT/2mn with fee de- 
noting the Boltzmann constant, T the temperature and 
ttih the mass of a hydrogen atom. 

2.3. Atmospheric opacities 

In this work, we adopt the dual-band approxima- 
tion for radiative transfer, namely that the black- 
body peaks of the stellar and exoplanetary emission 
are well separated in wavelength or frequency (see 
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Fig. 1. — The global-mean temperature-pressure profiles for Mod- 
els C, W and H. Models with 70 < 1 do not display temperature 
inversions, while models with 70 > 1 do. 



iHeng, Frierson feTP hillipps 20111 and references therein). 
One then requires two broadband opacities to completely 
describe the radiative transfer: a shortwave / optical opac- 
ity Kg to quantify the absorption of starlight with 
depth/pressure and a longwave/infrared opacity kl to 
quantify the absorption of thermal emission from the ex- 
oplan et. We neglect the effects of scattering (see iHengl 
2012). Varying Kg, while holding kl fixed, mimics the ef- 
fect of extra, shortwave absorbers in the atmosphere, of 
unspecified chemistry, which cause a temperature inver- 
sion when Kg > kl (Hubeny, Burrows, & Sudarksy 2003, 
Burrows et al. 2007b, Burrows & Orton 2010, Knutson 
et al. 2008). 

In the absence of empirical constraints, we assume the 
shortwave opacity to be constant, which means that it is 
related to the shortwave optical depth in a simple man- 
ner, 

k S P t So P 

T S = = "FT"* (6) 
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where Pq = 1 kbar is the reference pressure at the bottom 
of the simulation domain. The longwave optical depth 
is de scribed by two terms (jHeng. Frierson fc Ph illipps 

EoiiD, 
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TABLE 2 

Table of parameters, and their values, specific to each model 



Model 


f (W m- 2 ) 


n p (s- 1 ) 


a (AU) 


Unit (K) 


At (s) 


Simulation time (Earth days) 


Cim (J K- 1 m- 1 ) 


t- 1 is' 1 ) 


C 


2 X 10 4 


1.3 x 10~ 6 


0.28 


589, 501 


600 


4500 


10' 


2.1 x 10 _1 


CI 


5 x 10 4 


2.6 x 10~ 6 


0.17 


741, 630 


240 


3000 


10 7 


4.1 x 10 _1 


C2 


7 x 10 4 


3.4 x 10" 6 


0.15 


806, 686 


240 


3000 


10 7 


5.4 x 10 _1 


W 


2 x 10 5 


7.5 x 10~ 6 


0.09 


1048, 892 


180 


2000 


10 6 


1.2 


Wl 


5 x 10 5 


1.5 x 10~ 5 


0.06 


1317, 1121 


120 


1500 


10 6 


2.4 


W2 


7 x 10 5 


1.9 x 10~ 5 


0.05 


1433, 1219 


120 


1500 


10 6 


3.0 


H 


2 x 10 6 


4.2 x 10~ 5 


0.03 


1863, 1585 


120 


1500 


10° 


6.7 


HI 


5 x 10 6 


8.4 x 10~ 5 


0.02 


2343, 1994 


80 


1000 


10 8 


13.3 



Note: initial temperatures quoted are for 70 = 0.5 and 2.0, respectively. 



where the second term is an approximate way to account 
for the effects of collision-induced absorption, which in- 
troduces an enhancement factor of e at the bottom of 
the simulation domain. The longwave optical depth nor- 
malization is 7l = Ko-Po/5p, where kq is the longwave 
opacity normalization. We pick e = 1. Ignoring collision- 
induced absorption in the infrared is a good approxima- 
tion as long as this effect is only significant at atmo- 
spheric layers residing below the longwave photosphere 
(Hcng, Frierson & Phillipps 2011). 

In the absence of clouds/hazes, the ratio 70 = ks/^o 
determines if a temperature inversion exists in the at- 
mosphere (H ubeny, Burrows, fc Sudarksy 2003; see also 
discussion in lHeng et a,lj[201~2h . If 70 > 1, the shortwave 
photosphere sits above the longwave photosphere, and 
an inversion exists; if 70 < 1, the relative locations of 
the two photospheres are reversed, and there is no in- 
version. To study the effect of a temperature inversion 
on the day-night heat redistribution, we examine models 
with 70 = 0.5 and 2. We choose a typical value for the 
longwave opacity normalization: kq — 0.01 cm 2 g _1 such 
that the longwave/infrared photosphere resides at 
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^ = 0.1 bar. 
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The shortwave opacity follows from the chosen value of 
7o- 

The global-mean temperature-pressure profiles are 
showed in Fig. [TJ for the two chosen values of 70 = 0.5 
and 2, and the three representative flux strengths in- 
dicated with models C, W, H in Table [2] The effect of 
the temperature inversion is clearly evident: models with 
70 = 0.5 are hotter in the interior, and cooler in the out- 
erparts, while models with 70 = 2 display the opposite 
behaviour. We have set the internal heat flux of our hot 
Jupiter models to be zero. For completeness, we include 
in Appendix [XI the zonal-mean zonal wind, temperature, 
potential temperature and streamfunction profiles corre- 
sponding to these 6 models. 

2.4. Atmospheric dynamics 

The Princeton-GFDL FMS dynamical core 
solves the primitive equation s of meteorology (see 
iHeng. Menou fc Phillippsi 120111 and references therein). 
In the vertical, we choose the simulation domain to 
span 6 orders of magnitude in pressure, from P = 1 
mbar to 1 kbar; the latter value is chosen to roughly 
match the transition from the radiative to the convec- 
tive zone of a hot Jupiter. The vertical resolution is 
N v = 36, such that each pressure scale height is covered 
by about 3 grid points. The horizontal resolution is 
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Fig. 2. — Top: Ratio of night- to day-side hemispheric, infrared 
fluxes, for models with various degrees of irradiation. The night- 
day contrast decreases with Ti rr , indicating a weakening of the 
heat redistribution with increasing irradiating fluxes. Models with 
temperature inversions display a larger flux contrast, but this is a 
secondary effect compared to the irradiating flux. Bottom: Ratios 
between the advective and the radiative timescale, for the models 
of the top panel. Redistribution begins to break down when the 
advective timescale becomes much larger than the radiative time. 



192 x 96 (longitude versus latitude). In the absence of a 
physically-motivated scheme for more sophisticated ini- 
tial conditions, the simulations are started from a state 
of rest with a constant tem perature Ti^tH For e = 1, 
the i nitial temperature is (jHeng. Frierson fe Ph illipps 

1ml) 

16er S B \ 7o 



Ti 



init 



The specific value of T; n i t used to initiate each of the 
models is reported in Table [2] 

The areal heat capacity of the simulation bottom (Ci n t) 
is chosen to have a value such that the thermal iner- 
tia of gaseous hydrog en matches the radiative time con - 
stant at P = 1 kbar (|Heng. Frierson fc Phillipps! 120 111) . 
For Models C, W and H, it suffices to adopt Ci nt = 
10 7 , 10 6 a nd 10 6 J K" 1 m~ 2 , respecti v ely, a s it has been 
shown by He ng. Frierson fc Phillippsi (|2011| ) that order- 
of-magnitude variations in Cj nt produce small (< IK) 
variations in the temp erature-pressure profile. 

As already noted by IHeng. Menou fc Phillippsi (|2011[) . 

6 This particular initial condition is inconsequential for our simu- 
lations, since the isoth ermal initial condition in the code reproduces 
the radiative solution (Hcng, Frierson & Phillipps 2011). 
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Fig. 3. — Longitudinal separation between the location of the 
brightest spot on the exoplanet and the substellar point, i.e., the 
hotspot offset. Note that in the coldest models, the day-night con- 
trast is negligible, and the location of the brightest spot is subject 
to statistical fluctuations with time. When these conditions occur, 
we set the offset equal to 180°. 



numerical noise accumulates at the grid scale and has 
to be dissipated via a "hyperviscosity" in our spec- 
tral code. Since the hyperviscous term is neither an 
original nor a physical part of the primitive equa- 
tions, its value cannot be specified from first prin- 
ciples. As in IHeng. Menou & , Phillippsl (|2011[ ) and 
iHeng, Frierson fc Phillippsl (|2011l ). we pick the order-of- 
magnitude value of the hyperviscous timescale to be as 
large as possible in order for the simulations to collec- 
tively run to completion: t v — 2ir x lCP 6 /i! p . Decreasing 
the hyperviscous timescale — and hence increasing the 
rate of dissipation (t~ l ) — by one to two orders of magni- 
tude leaves the hotspot offset largely u nchanged (see Fig. 
19 of lHeng. Frierson fc P hillippsll 2 1 If) , but decreases the 
zonal wind speeds by ~ 10% (jHeng. Menou &: Phillippsl 
1201 If) . Our estimates for the zonal wind speeds are thus 
upper limits with regards to this numerical issue. 

Finally, we note that convective adjustment is applied 
when the dry lapse rate becomes super-adiabatic, i.e., 
the Schwarzschild criterion for convective stability is vi- 
ol ated. Our convective adjus t ment scheme is described 
in IHeng. Frierson fc Phillippsl (|2011f ). 

3. ATMOSPHERIC CIRCULATION, HEAT 
REDISTRIBUTION, AND HOTSPOT LOCATION 

In a strongly irradiated exoplanet, the atmospheric dy- 
namics and radiation are regulated by the interplay of 
zonal winds, which tend to redistribute/transport heat 
to cooler regions of the exoplanet, and local re-radiation 
of the absorbed energy, which favors the emission of ra- 
diation at the same location where it was absorbed. The 
relative efficiency of these two processes can be estimated 
by means of the associated timescales: 
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represents the fastest timescale on which zonal winds re- 
distribute temperature variations over the scale of the 
exoplanet, while the radiative timescale (|Goodv &: Yungj 
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is the characteristic timescale on which thermal energy 
is radiated away by a parcel of fluid. The latter is 
given by the ratio of the thermal heat content (cpT mSLX m 
where rh = P j g p is the column mass in hydrostatic equi- 
librium) to the re-emitted flux (csbT^^). Note that 
Jniax = T(raax.{v&}). In other words, both timescales 
are evaluated at the location where the zonal velocity 
attains its maximum value. Equation (|11[) is valid at op- 
tical depths close to or greater than unity; more precise 
versions of it carry an order of unity correction factor 
related to the efficiency of heat redistribution. 

Hot Jovian atmospheres possess radiative timescales 
which vary across many orders of magnitude. At the 
top of the atmosphere (~ 1 mbar), the atmosphere is 
radiative (t ra d -C i a dv); at the bottom (~ 1 kbar), it is 
advective (£ ra d 3> iadv)- The shortwave opacity controls 
the height/depth at which the bulk of the starlight is 
deposited. If the deposition occurs mostly near the top 
of the atmosphere (high ks), then the hotspot offset is 
small. Conversely, if it occurs near the bottom, then the 
offset is large. Thus, there is a degeneracy of the hotspot 
offset on irradiation and optical opacity. 

We denote the flux emerging from the longwave pho- 
tosphere by .Folr (where "OLR" denotes the "outgo- 
ing longwave radiation" , a term commonly used in the 
atmospheric/climate science community) . Denoting the 
latitu de by <I>, its latitudinal average is (jCowan fc Aeoll 
l2008h 

i r /2 

(^OLR) = - / J~ OLR COS 2 $ d$ . (12) 

7T J-tt/2 

With the longitude denoted by O, the total flux from 
each hemisphere of the exoplanet is defined as 



day, night 



(13) 



where {6 m j n , 6 max } = {t/2, 37r/2} for the day side. The 
night side is described by {© mm ,©max} = {0,7r/2} and 
{37r/2,27r}. The substellar point is set at O = n. In 
computing the day-night flux contrast, we have chosen to 
deal entirely with the latitudinally-averaged OLR or ID 
"brightness map" , rather than the thermal phase curve 
which is a convolution of the physical distribution of flux 
across longitude and its geo metric projection to the ob- 
server (jCowan fc Ago]ll2008r i. Additionally, we choose to 
plot the ratio of night-to-day fluxes, rather than the more 
familiar day-to-night ratio, since this is a quantity which 
ranges between and 1 (rather than between 1 and in- 
finity). However, our discussion will be mostly in the 
context of day-to-night flux ratios, which we generally 
refer to as the "day-night contrast" . 

Figure [2] shows the day-night contrast for the 16 ref- 
erence models that we have studied. The three coldest 
models C, CI and C2 display a near perfect redistribution 
of the heating, with comparable day and night fluxes. A 
slight contrast starts to appear for the W model (irradi- 
ated by JFo = 2 x 10 8 erg cm~ 2 s _1 ), and it then gradually 
increases with the strength of the irradiating flux. For 
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the same value of J-q, models with temperature inversions 
display a larger day-night flux contrast since the bulk of 
the starlight is deposited higher up in the atmosphere 
where the flow is dominated by radiative cooling rather 
than advection. The dependence of the flux contrast on 
opacity (and hence on the presence or not of tempera- 
ture inversions) is consistent with the results of Fortney 
et al. (2008) using ID models (see also Dobbs-Dixon & 
Lin 2008). 

The bottom panel of Fig. [2] shows the ratio between 
the advective and the radiative timescales as previously 
defined, for each of the models of the upper panel. Re- 
call that, since flow velocities vary with latitude, and so 
do the physical distances traveled by the flow from one 
side to the other of the exoplanet, the advective time 
is a function of $. To assign a unique value in Fig. [21 
we searched for the latitude at which the zonal flow at- 
tains its maximum speed, and computed the timescales 
at that latitude. For the coldest cases, the advective 
and the radiative times are within an order of magni- 
tude of each other. The fluctuations in £ a dv/irad arise 
from the fact that the difference between T max and the 
mean temperature of the flow becomes less pronounced 
as the strength of irradiation decreases. In other words, 
describing whether the zonal flow is predominantly ad- 
vective or radiative by a single number becomes a less 
robust exercise when the flow becomes more zonally sym- 
metric. Instead, the day-night flux contrast is a better 
indicator. 

As the irradiating flux becomes more intense, the ratio 
iadv/^rad rapidly increases. Our numerical results hence 
support the qualitative expectation that it is the rel- 
ative magnitude of these two timescales which mainly 
determines the extent of heat redis tribution, and conse- 
quent ly the day-night flux contrast ([Showman fc Guillotl 
2002). We note that, while our study has focused on 
spin-synchronized planets, the qualitative considerations 
derived above are expected to hold more generally. A 
change in the planet spin for the same T; rr would modify 
the flow dynamics (and hence both t rac i and t a dv); how- 
ever, the extent to which the flow is able to redistribute 
heat versus dissipate it locally, would still depend on the 
relative magnitude of these timescales, as long as the spin 
period is not much smaller than t rat \ (in which case the 
planet would be kept at a uniform temperature). 

The longitudinal offset ©offset between the hotspot lo- 
cation and the substellar point is displayed in Fig. [3[ 
For the three coldest models, the temperature is virtu- 
ally homogeneous on the entire exoplanet, but is sub- 
jected to statistical fluctuations; hence the location of 
the hottest spot is not well-defined. For these cases, we 
set 8 ffset = 180°. As the day- night contrast starts to 
deviate from unity, the hotspot location becomes well- 
defined. In the W model, the hotspot offset from the 
substellar point is rather large, ©offset ~ 130°-150°; this 
is indicative of efficient heat redistribution. As the irradi- 
ating flux increases, the hotspot location moves towards 
the substellar point. This can be easily understood in 
light of the changing ratio between the advective and 
the radiative timescales. As advection becomes less ef- 
fective in keeping up with the rapid cooling of the gas, 
the hotspot location moves closer to the point of maxi- 
mum irradiation, i.e., the substellar point. 

We note that the effects of opacity and mean molecu- 



lar weight were explored by iMenovJ (|2012aj ) using also an 
atmosph eric circulation code with dual-band radiative 
transfer (Rauschcr fe Menou I I2012T ) . albeit in the con- 
text of Neptune-like exoplanets. Our re sults are qualita- 
tively consistent with those reported in iMenoul (|2012af ). 
where it was found that variations in the shortwave opac- 
ity produce non-negligible differences in the hotspot off- 
set. Other studies of the role of atmospheric opacity 
in controlling the day / night contrast and dynamics have 
been performed by Lewis et al. (2010), tailored to the 
hot Neptune GJ436b. Their simulations particularly il- 
lustrated the important role of metallicity in regulating 
atmospheric circulation for a broad range of warm extra- 
solar planets. 

4. DISSIPATION WITHIN THE ATMOSPHERIC FLOW 

In this section, we use the 16 models previously de- 
scribed to investigate the importance of various forms 
of dissipation within the flow, and their relative trends 
with the strength of the irradiating flux. In particular, we 
first estimate the contribution to hydrodynamic dissipa- 
tion from gas which is either Kelvin-Helmholtz unstable 
and/or shocked, and then we study the additional contri- 
bution from Ohmic dissipation if the flow is magnetized. 

We remark that, since our simulations do not self- 
consistently account for the feedback effect of shocks, 
instabilities, and magnetic drag (all contributing to re- 
ducing wind speeds and hence the amount of dissipation), 
the estimates made in this section should be considered 
as upper limits for each particular combination of pa- 
rameters. However, the trends that we find, while likely 
to be less steep with the presence of feedback effects, are 
expected to be robust with respect to the dependence on 
irradiation. 

4.1. Hydrodynamic dissipation 

Our 3D circulation models, while not directly simulat- 
ing the onset and development of shocks and shear insta- 
bilities, allow us to establish whether the flow is subject 
to such instabilities and to shocks and to determine the 
size of the affected region. In the current study, we are 
primarily interested in exploring the effects of irradia- 
tion, while also examining the influence of (shortwave) 
opacity. 

4.1.1. Dissipation from Kelvin-Helmholtz-unstable gas 

The first form of hydrodynamic dissipation that we 
consider is due to the Kelvin-Helmholtz (KH) instabil- 
ity in the thermally-stratified shear flow. Whether such 
a flow becomes unstable can be estimated by measure- 
ments of the Richardson number 



TZi = N- 



where 



\ dz J 



N 



9 P de 7 
e T dz 



1/2 



(15) 



is the Brunt- Valsala frequency, 9? = T(P/ P )~ K is the 
potential temperature, and o"o = P/Po- The criterion 
for the flow to become Kelvin-Helmhol tz unstable is ex- 
pressed by the condition T Zi < 1/4 (|Kundu fc Cohen! 
l200ilLi fc Goodman ll2010ft . 
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Fig. 4. — Maximum Mach number as a function of depth within 
the atmosphere, for models with varying degrees of irradiation. 
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Fig. 5. — Top: Kinetic energy available within the flow, together 
with the amount available for dissipation through shear instabilities 
and shocks, for the three representative models C, W, H and no 
temperature inversions. Bottom: Same as above, but for the C, 
W, H models with temperature inversions. 



We computed the Richardson number at each point of 
the flow. For Model C, we find TZi > 1/4 everywhere, and 
hence shear instabilities are not expected at low levels 
of irradiation. In the warm and hot models, we find 
TZi < 1/4 only in the outermost pressure levels, down to a 
few millibars. Hence, we conclude from this analysis that 
the onset and consequent dissipation of shear instabilities 
is unlikely to play a major role in inflating hot Jupitcrs. 

4.1.2. Dissipation from shocked gas 

In order to estimate the presence and amount of dissi- 
pation due to shocked gas, we begin by keeping track of 
the Mach number at each location within the flow; since 
zonal velocites are the most significant ones, we define 
the Mach number as 

M(Q,*,r) = — . (16) 

C s 

Shocks are likely to occur for A4 > 1, barring special 
conditions which allow the pre- and post-shock region s 
to remain in causal contact (|Zel'dovich fc RaizerH '1966'). 
One may regard our simple estimates of the shock dis- 
sipation as upper limits since M. > 1 is a necessary but 
insufficient condition for the development of shocks. 

At each pressure level within the atmosphere, flow ve- 
locities change with both latitude and longitude. To in- 
vestigate how deep shocks can penetrate, in Fig. 2] we 
show the maximum Mach number as a function of atmo- 
spheric depth, for the three representative models C, W, 
and H, each with and without temperature inversions. 
In the coldest model, the flow is always subsonic. Wind 
speeds increase, as expected, with the strength of the ir- 
radiating flux and hence with the day-night temperature 
gradients. In the W model, the flow can become super- 
sonic (with Mach numbers of ~ 1.5 — 2) at pressure levels 
above a few tens of bars. In the H model, shock penetra- 
tion reaches a few bars, with Mach numbers ~ 2.5 — 3. 
For the same strength of irradiation, flows with temper- 
ature inversions reach a higher Mach number, due to the 
higher temperatures attained in the outermost layers. 

In order to make a crude estimate of the amount of en- 
ergy available to be dissipated in shocks, we first compute 
the total amount of kinetic energy available in the flow up 

to pressure P,E Mn (P) = (1/2) J*^ dVp[vl+v%). Note 
that the dominant term comes from the zonal component 
(i>e). Figure [5] shows -Ekin(-P) in the three representative 
models C, W and H. The cumulative amount from the 
whole flow varies by over three orders of magnitude be- 
tween the C and the H model, with a weaker dependence 
on the shortwave opacity. Generally, we find that the 
stronger the irradiation, the deeper the penetration of 
the zonal winds. 

Together with the local kinetic energy density, we 
track, at each point of the flow, both the Mach num- 
ber and the Richardson number. If either A4 > 1 or 
TZi < 1/4, we assume that a fraction of 0.5 of the ki- 
netic energy is dissipated/converted into heat. In strong 
shocks, this fraction is give n by 4(7 — + l) 2 , w here 
7 is the adiabatic index (|Draine fc McKei 119931) ; for 
7 = 7/3, this amounts to 0.48. For shear instabili- 
ties, Li & Goodman (2009) estimate the dissipated frac- 
tion of kinetic energy to be « 0.5-1. Hence we take 

E dis (P) = (l/2)/dis/; ( p ) ) dVp{v%+v%), where / dis = 0.5 
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Fig. 6. — Dissipated Ohmic power as a function of depth, for the three representative models C, W, H, both with and without temperature 
inversions. A magnetic field strength of 2 G at the poles has been assumed. Note that the penetration depth increases with the strength 
of the irradiating flux; for the same T; rr , it is larger in models without temperature inversions. 



if A4 > 1 or TZi < 1/4 and /di S = otherwise. We remark 
once again that this is a rather crude estimate, motivated 
by the fact that current global 3D circulation models are 
not yet able to follow the development, growth, and dis- 
sipation of instabilities and shocks. 

The cumulative dissipated energy, Edi S {P), is shown 
in Fig. [5] for the C, W and H models. In the C model, 
Edis(P) = since, as discussed above, the criterion for 
the onset of the KH instability is never satisfied, and the 
velocities are always subsonic. In the W and H models, 
on the other hand, a fraction of about 10-20% of the ki- 
netic energy is available for dissipation, with the fraction 
being generally larger in models with temperature inver- 
sions, due to the stronger velocities (and their gradients) 
in the upper atmospheric layers, which are the ones more 
prone to becoming either KH unstable or shocked. For 
a given value of 70, dissipation occurs down to deeper 
levels as the strength of the irradiating flux increases, re- 
flecting the penetration depth of the shocks (cf. Fig. |4|). 
However, even in the hottest models, shock dissipation 
never penetrates deeper than a few bars. 

4.2. Ohmic dissipation 

Previous work (Batygin et al. 2010; Perna et al. 
2010b) has shown that, for Jupiter-like magnetic field 
strengths, Ohmic dissipation can be a significant source 



of energy for the exoplanet. In particular, using a simu- 
lation specialized to the physical parameters and con- 
ditions of HD 209458b, Perna et al. (2010b) showed 
that Ohmic dissipation can penetrate down to levels deep 
enough to influence the radius evolution of the exoplanet. 
Since irradiation is the driver of the flow dynamics, cur- 
rents (and hence Ohmic dissipation) are expected to de- 
pend sensitively on its strength. In the following, we 
investigate this dependence using our 16 baseline mod- 
els. We assume the £?-field to be an axisymmetric dipolc 
where the rotational and dipolar axes are aligned. The 
dependence of Ohmic dissipation on the magnitude of the 
£?-field has been explored (inclusive of an approximate 
feedback from magnetic drag), by Perna et al. (2010b). 

To compute Ohmic dissipation within the context of 
our simulations, we follow the formalism of Liu et al. 
(2008) and Perna et al. (2010b). Under the assump- 
tion that zonal winds are stronger than the meridional 
motions, the only relevant component of the induction 
equation is the toroidal one. Under steady-state condi- 
tions, it was showed by Liu et al. (2008) that the result- 
ing dominant component of the current induced by the 
zonal flow is the meridional one, which we write below 
in terms of the latitude $ and the longitude 9, and also 
incorporating the expression for an aligned dipole for the 
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Here, rj = c 2 /(Aira e ) is the electrical resistivity, and the 
electrical conductivity <j e is computed using the same 
formalism as in Perna et al. (2010a, b). The familiar ra- 
dial coordinate, denoted by r, is equal to z + R p where 
z is the vertical height measured from the bottom of the 
simulation domain. The top of the simulation domain is 
located at r max = max(z + R p ). We choose to compute 
the current density as a function of the variable z, which 



is defined such that z — when r 



i.e., it is the 



vertical height measured from the top of the model at- 
mosphere. Again, we is the zonal velocity. The quantity 
Bq is the normalization term in the dipole magnetic held, 
such that it is B at the equator and 2B at the poles. 
For our calculations, we choose B = 1 G. Note that in 
the above equation we have neglected the contribution 
from an unknown boundary current. 

The total Ohmic power dissipated between the bottom 
shell at pressure Pq and a shell of material located at 
pressure P, is then readily computed as 



Pj(P) = / dr'r 

Jr(P) 



2tv 



dO 



d® cos $ 



er e (r',©,$) ' 
(18) 

Fig.[6]shows the cumulative Ohmic power as a function 
of the atmospheric depth, for the three representative 
models C, W and H, both with and without tempera- 
ture inversions. The magnitude of Pj has a very strong 
dependence on the irradiating flux, varying at maximum 
by about 16-17 orders of magnitude between the C and 
the H models. While the increasing zonal wind speeds 
with T; rr do contribute to this effect, the dominant factor 
is the exponential dependence of the electrical conductiv- 
ity on temperature. 

For the same Ti rr , the maximum Ohmic power is larger 
for 70 = 2 (with a temperature inversion) than it is for 
7o = 0.5 (no inversion), while the penetration depth on 
the other hand is larger for 70 = 0.5. This behaviour 
can be understood as follows. For 70 = 2, the outer 
layers of the atmosphere are hotter than the inner lay- 
ers, and hence the conductivity is much larger in the 
outer atmospheric regions. Since velocity gradients are 
also stronger in these regions, the net result is the gen- 
eration of very large peripheral currents and hence large 
dissipation, which however saturates, even in the hottest 
models, close to the photosphere. For 7 = 0.5 on the 
other hand, the outer layers are colder, and the low con- 
ductivity results in weak currents. Deeper in, the tem- 
perature rises, and so does the conductivity, resulting in 
an increase of the Ohmic power at deeper atmospheric 
levels. However, the maximum dissipated power satu- 
rates at smaller values than than for the case of 70 = 2, 
since for 70 = 0.5 the region with higher conductivity 
corresponds to zonal flows with weaker gradients. 

A comparison between the total power (i.e., integrated 
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Fig. 7. — Comparison between hydrodynamic and Ohmic dis- 
sipation, for all of the models presented in our study. For the 
hydrodynamic dissipation, we assume that the energy is dissipated 
over a rotational timescale. The Ohmic power is computed for a 
dipolar _B-ficld with strength of 2 G at the pole. Without feedback 
from magnetic drag, the Ohmic power scales as B?. The dotted 
lines labelled and "10%" correspond to 1% and 10% of the 

stellar irradiation, respectively. 

over the entire atmosphere) generated through shocks 
and shear instabilities, and through Ohmic dissipation, 
is shown in Fig. [7J Since our simulations do not explic- 
itly follow shocks and the development of shear insta- 
bilities, they cannot make predictions for the timescales 
over which dissipation associated with these phenomena 
occur. To compute the hydrodynamic power in Fig. 
we assumed that -Edis is dissipated over a rotational 
timescale (which is comparable in magnitude to the ad- 
vection timescale). The 2D simulations by Li & Good- 
man (2010) show that, once shocks are formed, they slow 
down over a very short time (shorter than the advection 
time). As long as the dissipation timescale is shorter 
than 27r/f2 p , the power shown in Fig. [7] can be consid- 
ered as a good proxy for its corresponding time-averaged 
quantity over a rotational cycl^H- 

Our simulations show that the Ohmic power is a much 
steeper function of irradiation than the hydrodynamic 
power, which we find to be dominated by shocks. For 
the latter, only velocities matter, while the dissipated 
Ohmic power depends on velocity gradients. But most 
importantly, as the temperature increases, so does the 
ionization fraction, and hence the conductivity, which 
results in stronger currents. However, we remark, once 
again, that the slope at high T irr would be shallower 
if magnetic drag were incl uded se l f-consis tently in the 
simulations. Specifically, iMenoul (|2012bl ) has shown, 
via simple scaling laws, that magnetic drag rapidly re- 
tards the atmospheric winds for strong stellar irradia- 
tion, and that, for a given magnetic field strength, there 
is a value of the irradiation/equilibrium temperature for 
which Ohmic dissipation peaks (and then declines as ir- 
radiation becomes more intense). Our results for the 

7 Lacking 3D evolutionary planet models which can assess the 
effect of anisotropic heat injection in the planet interior on the 
thermal evolution of the planet, for an order-of-magnitude com- 
parison with the ID evolutionary models of Guillot & Showman 
(2002), we consider the averaged dissipated hydrodynamic power 
over a rotational cycle. 
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dissipated power should therefore be considered as up- 
per limits. 

Finally, we note that iHeng I ((201 2T ) provides a compan- 
ion study to our paper by using semi-analytical models 
to elucidate the interplay between atmospheric scatter- 
ing, absorption and Ohmic dissipation with relevance to 
the possible presence of atmospheric clouds or hazes. 

4.3. Implications for the problem of radius inflation 

As discussed in §1, the anomalously large radii of a 
fraction of hot Jupiters remains a persistent problem 
for researchers. An "extra" source of internal energy is 
needed for the radii to be larger than predicted by stan- 
dard evolutionary models. Guillot & Showman (2002) 
argue that, if about 1% of the stellar insolation flux were 
deposited at pressures of tens of bars deep in the atmo- 
sphere, cooling would slow down sufficiently to explain 
the inflated radii of hot Jupiters. The evolutionary model 
presented by Guillot & Showman (2002) was tailored to 
the specific case of the exoplanet HD 209458b, with a 
radius of R = 1.35i?j is the Jupiter radius), and 
a convective-radiative boundary located at pressures of 
about 1 kbar (at its current age). Note that, for exoplan- 
ets with comparable mass, the boundary layer between 
the convective/radiative zones is located at lower pres- 
sures as the planet is younger and larger (e.g., Guillot 
& Showman 2002, Burrows et al. 2003). The amount 
of internal heat needed to inflate the radius by a certain 
AR is also a strong function of the planet mass (Miller, 
Fortney & Jackson 2009). The smaller the mass, the eas- 
ier it is to inflate the planet. For example, Miller et al. 
(2009) show that a 20% increase in radius (with respect 
to i?j) can be achieved with a mere 10 24 erg s" 1 in a 
planet of mass 0.1 Mj, while the same radius increment 
would require about 10 28 erg s _1 for a 10 Mj planet. 

While several mechanisms could transport and dissi- 
pate energy into the interior, such as eccentricity damp- 
ing (Bodenheimer et al. 2001), or forced turbulent mix- 
ing in the radiative layer (Youdin & Mitchell 2010), here 
we have investigated energy deposition resulting from the 
development of shear instabilities and shocks, and from 
Ohmic dissipation. For all of these phenomena, our nu- 
merical simulations, while not modeling the dissipation 
itself and its feedback on the flow, do allow us to esti- 
mate their importance from an energetic point of view, 
as a function of the strength of the irradiating flux and 
of the atmospheric opacity. 

We have found that shear instabilities can develop 
only very superficially, even in the hottest models, and 
hence they are unlikely to meaningfully influence the 
energy budget of the exo planet. If the source of w eak 
turbulence postulated by lYoudin fc Mitchell I (|2010f ) re- 
sults from shear instabilities, then their proposed mecha- 
nism of forced turbulent mixing does not operate deeply 
enough. To further test the relevance of their mechanism 
in global circulation models requires the specification of 
an instability and its criterion for being invoked. 

Shocks penetrate deeper, up to a few bars in the hottest 
models. This is a depth which could be interesting for 
influencing the thermal evolution of planets, and hence 
the degree by which their radius is able to shrink during 
the evolution. For the particular case of HD 209458b, 
the evolutionary models of Guillot & Showman (2002) 
showed that, even at a pressure as shallow as 5 bars, 



dissipation of 10% of the absorbed stellar flux (corre- 
sponding to E = 2.4 x 10 28 erg s _1 for this planet) is 
sufficient to retard cooling so that the planet contracts 
to its current size on a timescale ~ 5 Gyr, comparable to 
its age. Inspection of the trend for the dissipated hydro- 
dynamic power versus irradiating flux in Fig. shows 
that, for the strength of Jo corresponding to the case 
of HD 209458b, the hydrodynamic power is on the or- 
der of a few xlO 28 erg s _1 . Therefore, since the dissi- 
pated hydrodynamic power that our models predict falls 
in an interesting order of magnitude range (at least for 
the most irradiated objects), and its penetration depth 
also reaches levels close to the interesting ones identified 
by evolutionary models^, we suggest that shock dissipa- 
tion could indeed play a significant role in the thermal 
evolution of the planets, and hence influence their radii, 
although especially so in the most irradiated planets. 

Our investigation of Ohmic dissipation has showed 
that, for a -B-field strength on the order of a Gauss, the 
total dissipated Ohmic power is generally lower than the 
hydrodynamic power except for the most irradiated hot 
Jupiters. However, it penetrates deeper, down to pres- 
sures of several tens of bars, if the atmosphere has no 
temperature inversion. In the evolutionary scenario of 
Guillot & Showman (2002), a deeper penetration depth 
can be "traded" for a smaller power. For example, in the 
case of HD 209458b, their evolutionary models could re- 
produce the exoplanetary radius at its current age both 
with 10% of the irradiating flux dissipated at 5 bars, or 
with 1% (i.e., 2.4 x 10 27 erg s~ x ) dissipated at 21 bars. 
For a S-field strength of a few Gauss, this power is avail- 
able in the most irradiated objects. Clearly, the impor- 
tance of Ohmic dissipation for the thermal evolution of 
planets also depends strongly on the magnitude of the 
B- field, since the Ohmic power scales aqj Bq, and, for 
fixed planet conditions, the penetration depth also in- 
creases with larger Bq (Perna et al. 2010b). However, 
for a given i?-field, since both the penetration depth and 
the power dissipated at a fixed pressure are rapidly in- 
creasing functions of the strength of the irradiating flux, 
the influence of Ohmic dissipation on radius inflation is 
expected to depend strongly on stellar irradiation. 

5. QUALITATIVE COMPARISON WITH OBSERVATIONS 

Observations of hot Jupiters have allowed to character- 
ize, albeit with varying statistical level, the atmospheric 
properties that we have discussed so far, and in partic- 
ular, the efficiency of heat redistribution, the hotspot 
offset and the exoplanetary radius. Here we do not at- 
tempt to make any detailed, quantitative comparison be- 
tween our models and the data; the former have been 
constructed for "typical" values of the exoplanetary pa- 
rameters, and hence they are not generally applicable to 
a specific object. For the same value of the irradiating 

8 It is also important to note that vertical advection can 
help transport the dissipated energy to deeper levels. In our 
models, the vertical advection timescale is ~ 10 5 s, compara- 
ble to the advection timescale; for an HD 209458b-like planet, 
Heng. Fricrson & Phillipps (2011) showed that vertical mixing ex- 
tends down to about 10 bars. 

9 As discussed in §4.2, this scaling would be shallower if feedback 
effects from magnetic drag were accounted for, but it would still 
remain an increasing function of Bo for a wide range of irradiating 
fluxes. 
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flux, a scatter is expected due to the different intrinsic 
properties of the exoplanet, of which here we have only 
examined the d ependence on the shortwave opacity for 
two values. (See lHengll2012l for a broader exploration of 
70 values, and also the influence of shortwave scattering, 
using semi-analytical models.) The data on the other 
hand, still suffer from measurement uncertainties which 
would make quantitative comparisons with models not 
fully reliable yet. However, the qualitative trends in both 
the data and the model are fairly robust, and these are 
the ones that we are interested in reproducing at this 
stage. 

5.1. Data summary 

Figure [8] summarizes the main observational con- 
straints on the radius inflation and atmospheric circu- 
lation of hot gas giant exoplanets. The exoplanetary ra- 
dius measured for transiting gas giants (M > O.SMj) 
is shown as a function of the irradiation temperature. 
Different symbols indicate more/less massive exoplanets, 
exoplanets with or without a temperature inversion mea- 
sured near the infrared photosphere, and exoplanets with 
or without efficient day-night temperature redistribution. 
The presence of a temperature inversion is inferred from 
the broadband infrared emission spectra measured dur- 
ing secondary eclipse with the Spitzer satellite. The effi- 
ciency of the day-night temperature redistribution is in- 
ferred from the comparison of the effective temperature 
in the Spitzer passbands with the equilibrium tempera- 
ture. If the measured temperature is much larger than 
the equilibrium temperature, the redistribution of heat 
to the night side is likely to be lo\4 TO l. In one case, HD 
189733b, the heat redistribution is measured more di- 
rectly, from the infrared phase curve of the exoplanet at 
8 and 24 microns (Knutson et al. 2007, 2009). 

The exoplanet radii and irradiation data have 
been taken from the edited compilation at 
www. inscience . ch/transits (see references therein for 
primary sources). Information on the absence/presence 
of temperature inversions is adopted from Knutson et 
al. (2011, see references therein), while information 
on redistribution is taken from the compilation of 
Harrington (2011), updated with more recent data from 
Spitzer for HD80606, XO-3, HAT-P-7 and WASP-12; we 
also include inferences on XO-4, HAT-P-6 and HAT-P-8 
(Todorov et al. 2011). It is important to note that 
inferences on redistribution (presented in the lower 
panel of Figure 8) are based on data that stretches 
the capacity of the Spitzer instruments to their limits, 
and that the uncertainties associated with instrumental 
effects are still high. In at least two cases, Ups And b 
and HD 149026b, further measurements have shown 
the initial estimates for the surface temperature and 
redistribution efficiency to be completely incorrect. 
These two objects are not in our plot, Ups And b 
because it is not a transiting planet and therefore its 
radius is not measured, and HD 149026b because it is 
not in synchronous rotation and is a core-dominated 

10 This inference is tentative, since the albedo and the redistri- 
bution of flux in frequency by molecular spectral features will also 
play a role. The indications for specific objects may turn out to 
be incorrect, but hopefully this will simply introduce some random 
noise and not entirely mask the existing correlations between the 
efficiency of redistribution and other parameters. 



planet rather than a hot Jupiter, and therefore it does 
not fit our model assumptions. 

Information on the angular offset of the hotspot is only 
sparsely available. Spitzer observations of HD 189733b 
show that the hottest region of the thermal phase curve 
is advected to the east of the substellar point (Knutson 
et al 2007) by 16 ± 6 degrees; Knutson et al. (2009) 
inferred the ID brightness maps at 8 and 24 microns 
and found the offsets to be 30 ± 4 and 23 ± 7 degrees, 
respectively. Measurements for a couple of other objects 
are less certain. Crossfield at al. (2010) report, for Ups 
And b, which has a good redistribution, an offset of the 
ID brightness map to the east of either 57 ± 21 degrees 
or 84.5 ± 2.3 degrees, using two different data sets for the 
analysis. Cowan et al. (2012), report, for WASP-12b, an 
offset to the east of 12 ± 6 degrees from the Spitzer phase 
curve at 4.5 microns, and an offset of 53 ± 7 degrees at 
3.6 microns. 

5.2. Observational vs. theoretical trends 

Figure [8] illustrates the strong correlation between irra- 
diation and radius inflation. This indication that injec- 
tion of a fraction of the irradiation into the exoplanet as 
internal entropy is responsible for the radius inflation (as 
first suggested by Showman & Guillot 2002) is reinforced 
by the observed correlation between radius inflation and 
the mass of the exoplanet (Fortney et al. 2011), and by 
the absence of radius inflation for cooler Jupiter-sized ex- 
oplanet candidates from the Kepler mission (Demory & 
Seager 2011). Hence, the observations strongly support 
explanations of radius inflation as a consequence of high 
stellar irradiation. The observed slope of the dependence, 
though, seems even steeper than that predicted by the 
injection of a constant fraction of the stellar flux as inter- 
nal entropy as proposed by Showman & Guillot (2002), 
suggesting that the efficiency of the coupling mechanism 
increases with temperature. There is no significant cor- 
relation between radius inflation and orbital eccentricity 
(Pont 2009), nor with temperature inversion. 

For the dissipation mechanisms that we have explored 
in this work, both of hydrodynamic and of magnetohy- 
drodynamic origin, our simulations show a strong corre- 
lation with the strength of the irradiating flux. In par- 
ticular, for a B-ficld on the order of a Gauss, we found 
that, for Ti rr < 2000 K, the dissipated Ohmic power is too 
small (and its penetration depth too shallow), to affect 
the radius evolution in any appreciable way. However, as 
Tj rr approaches ~ 2000 K, the Ohmic power grows notice- 
ably, and its penetration depth, reaching down to several 
tens of bars, can then affect the evolution of the exoplan- 
etary radius (Guillot & Showman 2002). Interestingly, as 
Fig.[8]shows, radius inflation starts to become important 
as the irradiation temperature approaches the 2000 K 
range. We also note that variations of the magnetic 
field from exoplanet to exoplanet can result in a sub- 
stantial scatter in the radius evolution, since the Ohmic 
power scales as Bq , and this rapid growth with magnetic 
field strength dominates over the reduction to the power 
caused by the feedback effect of magnetic drag, at least 
for an interesting range of B-field strengths (Perna et 
al. 2010b). For a Jupiter-like magnetic field (~ 14 G at 
the pole), the power produced by Ohmic dissipation can 
already affect planets with T irr ~ 1700-1800 K. The at- 
mospheric (shortwave) opacity further contributes to the 
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Fig. 8. — Upper panel: Exoplanetary radius as a function 
of irradiation temperature for transiting gas giants. Exoplan- 
cts in the range 0.3 < M < l.5Mj are shown as closed sym- 
bols, exoplanets heavier than 1.5 Mj as open symbols (data from 
www.inscience.ch/transits). Dotted lines show the temperature 
range where no significant radius anomaly is observed (flat line), 
and least-square linear fits to the low-mass data and high-mass 
data separately. Middle panel: Same data as in the upper panel, 
but with upside-down red triangles indicating inverted temperature 
profiles, and blue triangles normal profiles. Bottom panel: Same 
data as in the other two panels, but with blue diamonds indicating 
probable efficient redistribution, and red circles probable strong 
day-night contrast. 



scatter, since the penetration depth of Ohmic dissipation 
is a non-negligible function of this variable. 

In summary, and at a qualitative level, our simula- 
tions predict that both hydrodynamic dissipation (dom- 
inated by shocks) and Ohmic dissipation (for magnetic 
field strengths at least as large as several Gauss), can 
play an important role in altering the thermal evolution 
of hot Jupiters, slowing down cooling, and favoring the 
presence of planets of larger dimensions than predicted 
by evolutionary models without "extra" energy injection. 
A significant dissipated power (both hydrodynamic and 
Ohmic if Bq > 7 — 8 G), on the order of several per- 
cents of the stellar insolation, is predicted for irradiation 
temperatures T irr > 1700 - 1800 K. For both forms of 
dissipation, the stronger the irradiating flux, the larger 
the total dissipated power and the deeper the penetration 
depth, resulting in an increasingly larger inflation as the 
flux becomes more intense. At large insolation, even if 
the direct dissipation depth is rather shallow (as it is for 
shocks and for currents in highly inverted temperature 
profiles), vertical mixing can help to carry down heat to 
levels of tens of bars, hence maintaining the energy in- 
jection in the relevant regions for affecting the thermal 
evolution of the planet. 

A quantitative prediction of the slope of the R p (Ti rl ) 
curve would need a coupling between a specific dissipa- 
tion mechanism and an evolutionary model for the exo- 
planet, which is beyond the scope of the present study. 
But even so, as discussed above, a large scatter would be 
expected due to exoplanets with different intrinsic prop- 
erties (opacity, gravity, magnetic field strength, etc). 

Although the statistical information is still rather lim- 
ited, inspection of Fig. [8] further shows that the effi- 
ciency of day-night heat redistribution seems to dimin- 
ish rapidly for T; rr above 2200 K, with the hottest gas 
giants showing day-side temperatures in the infrared 
much higher than the equilibrium temperature, an ef- 
fect pointed out by Harrington (2011). Interestingly, our 
simulations (cf. Fig. [2]) also show that, from irradia- 
tion temperatures around 2200-2500 K, the atmospheric 
flow pattern begins to evolve from a very good redis- 
tribution to a regime of no-redistribution. Exoplanets 
with temperature inversions evolve more quickly towards 
a no-redistribution regime by displaying, for the same ir- 
radiating flux, a larger day-night flux contrast. This is 
again supported by observations, as well as by other work 
(Fortney et al. 2008). It may also be worthwhile to note 
that the drop in redistribution efficiency appears to occur 
around the temperature at which radius inflation satu- 
rates, although this is only a tentative indication at this 
point due to the low number of known cases, the large 
scatter at high temperatures, and the assumptions in- 
volved in inferring redistribution from infrared day-side 
fluxes in the Spitzer passbands. From a theoretical point 
of view, a saturation of radius growth at very high values 
of T irr is expected due to feedback effects (e.g., shock dis- 
sipation and magnetic drag, both contributing to reduce 
the speeds of the fastest flows). 

For the hotspot offset with respect to the substellar 
point, there is insufficient data statistics yet to compare 
theoretical and observational trends. However, we do 
note that, as the observations show that the offset can be 
both rather small or very large, our theoretical modeling 
also shows a wide range of values, and further predicts a 
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correlation between offset width and degree of redistribu- 
tion. Although our modeling here has been for "typical" 
exoplanets and is not tailored to the properties of specific 
objects, it is however tempting to point out that, with 
the irradiating flux level of HD 189733b, our predicted 
offset is in the 20°-40° range (cf. Fig. [3]), consistent with 
the range of values (16°-34°) obtained from the inferred 
ID brightness maps, at 8 and 24 microns, by Knutson et 
al. (2009). In the case of Wasp-12b, the extremely high 
value of the irradiation temperature, 2j rr ~ 3500 K, is 
above the temperature range studied in our simulations. 
However, inspection of the offset trend at high T; rr in 
Fig. [3] would suggest a value around a couple of tens 
of degrees at those high temperatures. The measured 
offset of 16 ± 4 degrees at 4.5 microns appears compat- 
ible with this estimate. On the other hand, the much 
larger offset of 53 ± 7 degrees at 3.6 microns does not. 
The difficulty in reconciling this large phase offset with 
the other properties of the object was already discussed 
by Cowan et al. (2012), though they also noted that 
the measurement could be affected by highly-correlated 
residuals near the purported peak. If real, such a mea- 
surement might indicate a rather unusual opacity, and a 
detailed theoretical modeling would require a non-grey 
radiative transfer treatment (Fortney et al. 2006, Show- 
man et al. 2009, Burrows et al. 2010). Similarly, for 
the case of Ups And b, the unusually large offset values 
reported by Crossfield et al. (2010), albeit with some- 
what large error bars, cannot be simply accounted for 
with the opacity values adopted here. However we note 
that, even within the simple framework of the dual-band 
model, larger offsets can be produced for smaller values 
of both the optical and infrared opacities. 

6. SUMMARY 

Irradiation by the parent star is the main driver of 
the radiative and dynamical properties of hot Jupiters. 
Opacity in the exoplanetary atmosphere further plays an 
important role in determining the response of the atmo- 
spheric flow to a certain amount of irradiating flux. 

In this paper, by means of 3D atmospheric circulation 
models with dual-band radiative transfer, we have ex- 
plored the role of irradiation and opacity in determining 
some of the main observational properties of hot gaseous 
exoplanets, with particular emphasis on heat redistribu- 
tion and on hydrodynamic and magnetohydrodynamic 
dissipation. Although our simulations are still idealized 
in several respects (simplified radiative transfer, constant 
opacities, lack of feedback effects due to dissipation of in- 
stabilities and shocks and to magnetic drag), they still 
allow us to capture some broad trends exhibited by the 
observations. 

Our results are summarized below: 

• For irradiation temperatures T- m < 2200 — 2400 K, 
we find heat redistribution to be very efficient. 
Near-perfect redistribution (i.e., Fd ay /-Flight ~ 1) 
only occurs up to T irr < 1500 K, but in the range 
1500 K < T irr < 2400 K the day-night flux contrast 
is still relatively small (< 2). As the stellar irra- 
diation becomes more intense, redistribution be- 
gins to fail. Observations, albeit limited by the 
low statistics, hint at the redistribution breaking 
down around T- m > 2200 — 2400 K, in broad agree- 



ment with our theoretical results. Our simulations 
support the simple intuitive result that it is the 
interplay between advection and radiative cooling 
which determines the effectiveness of heat redis- 
tribution. Redistribution begins to falter as the 
advective timescale becomes much longer than the 
radiative timescale. 

• For the same strength of irradiation, exoplanets 
with a temperature inversion (here parameterized 
by means of a larger shortwave opacity) display a 
higher day-night contrast than exoplanets with no 
temperature inversion, since starlight is deposited 
higher up in the layers of the atmosphere domi- 
nated by radiative cooling (rather than advection) . 

• The offset of the hottest region (hotspot) from the 
substellar point can exhibit a wide range of values; 
it is generally very large (even entering the night 
side) for exoplanets with very efficient redistribu- 
tion (low T; rr ) and it moves closer to the substel- 
lar point with increasing strength of the irradiating 
flux. For a fixed value of the irradiating flux, opac- 
ity further plays a role in determining the offset. 
The smaller the shortwave opacity, the larger the 
offset for the same T; rr . 

• A measurement of the Richardson number through- 
out the 3D flow showed that the Kelvin-Helmholtz 
instability is unlikely to play a major role in the 
thermal evolution of the planet. We found that the 
the criterion for the instability onset, TZi < 1/4, is 
satisfied only in the uppermost layers of the flow, 
down to several millibars (if at all). 

• Shocks, associated with the presence of supersonic 
flows (up to Mach numbers of ~ 2.5 — 3 in the 
hottest models), can dissipate up to 10-20% of the 
available kinetic energy. The total dissipated power 
over an advective timescale reaches several percent 
of the insolation energy power for irradiation tem- 
peratures T irr > 1700 - 1800 K, and it then in- 
creases with Tj rr . Hence, for the most irradiated 
objects, the magnitude of the hydrodynamic power 
is such that it could noticeably influence the ther- 
mal evolution of planets. However, the direct pen- 
etration depth of the dissipation does not extend 
beyond a few bars even in the most irradiated hot 
Jupiters. Some amount of vertical advection would 
then be necessary to bring the heat to deeper lev- 
els. In the most irradiated objects, our simulations 
predict vertical mixing to penetrate down to tens 
of bars (see Appendix EJ, which could make hy- 
drodynamic dissipation an interesting contributor 
to the thermal evolution of the exo- Jupiters, and 
hence affect the evolution of their radii. 

• For a Jupiter-like B-field strength, Ohmic dissi- 
pation is negligible up to irradiation temperatures 
Tin- ^ 1700— 1800 K, as far as the power needed to 
influence the thermal evolution of the hot Jupiters. 
As the insolating flux becomes stronger, the mag- 
nitude of the dissipation increases, and so does its 
penetration depth. The latter effect is however 
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much more pronounced in models without temper- 
ature inversions, which have a hotter interior. Dis- 
sipation can occur down to levels of several tens of 
bars, hence potentially affecting the radius evolu- 
tion of the exoplanet. If Ohmic dissipation does 
play a major role in radius inflation, a strong cor- 
relation with the strength of the irradiating flux 
is naturally predicted. However, a large scatter is 
also expected, due to the strong dependence of the 
Ohmic power on the magnetic field strength, as well 
as on the atmospheric opacities. 
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APPENDIX 

ZONAL WIND, TEMPERATURE, POTENTIAL TEMPERATURE AND STREAMFUNCTION PROFILES 

For completeness, we include the temporally- and zonally-averaged zonal wind (Figure |9|), temperature (Figure [TT)|) . 
potential temperature (Figure lll|) and Eulerian mean streamfunction (Figure I12j) profiles as functions of latitude 
and pressure/height. The con tours of constant potential t empe rature are equivalent to contours of constant specific 
entropy. As already noted in iHeng. Frierson fc Phillippsl (|2011l) , computing the Eulerian mean streamfunction only 
on the dayside hemisphere strictly violates the definition of a streamfunction, but it does illustrate the nature of the 
global circulation pattern and the depth to which it penetrates. 
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Model C, 7 =0.5 




Model C, 7 =2.0 




Fig. 9. — Temporally- and zonally-averaged zonal wind profiles for Models C, W and H and for 70 =0.5 and 2. Contours are in units of 
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Model C, 7 =0.5 




Model C, 7 =2.0 




Fig. 10. — Same as Figure [9] but for the temporally- and zonally-averaged temperature profiles. Contours are in units of K. 
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Model C, 7 =0.5 




Model C, 7 =2.0 




Fig. 11. — Same as Figure[9]but for the temporally- and zonally- averaged potential temperature profiles. Contours are in units of K. 
Lines of constant potential temperature are equivalent to lines of constant entropy (isentropes). 
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Model C, 7 =0.5 (day side) 




Model C, 7 =2.0 (day side) 
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Model H, 7 =0.5 (day side) 
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Model H, y =2.0 (day side) 




Fig. 12. — Same as Figure[9]but for the temporally- and zonally-averaged Eulcrian streamfunction profiles only on the dayside hemisphere. 
Contours are in units of 10 kg s" 1 . Negative and positive values indicate clockwise and anti-clockwise circulation, respectively. 



